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THE DERIVATIVE OF AN IDENTITY 


WITH THE ANALYSIS OF FUNCTIONS OF VECTOR STATES 


WARNING: Do not try to create a Geometric picture while studying this material. The depth of the Mathematics is constructed 
from a rigorous application of the Hermitian CONJUGATE VARIABLE, and NOT functions of two variables. 


PREREQUISITE 

The Mathematics in this text is of High Order Source(Advanced). Its beyond High School, and Tertiary. An Adequate 
background on the concepts of the following is desired: 

Characteristic Polynomial Function 

Differentiation Rules 

Sequences and Series 

Partial Derivatives 

Linear Approximation 

The Average Value of a Function 

Standard Equation of the Circle with Radius (n) 

Vectors 

Complex Numbers 
Tensor Calculus 
Pauli Matrices 
Integration 

Programming Language 


PREFACE 


It is a bad habit for Physicists to take their most successful abstractions to 
be real properties of the world - Nathaniel Mermin 

To the Mathematician, Computer Programmer, Electrical Engineer and Physicist, 

Understanding the concept of an Absolute Value of State Space, is the most gratifying experience, and 
an achievement that takes lots of intense study and patience. 

The title of my previous book [Rational Commensurability and Why Division by Zero Is Not Possible] 
pointed on the properties of the number 0; while this text now REFLECTS the number 1. I must admit 
the structure of both materials are complete that I could not add or subtract any formalism to their 
architecture. 
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The depth of Mathematical maturity in visualizing its core, is the motivation for the title, The Derivative of 
An Identity. 

I’m most humbled to be able to clarify a well and thoroughly established principle. 

INTRODUCTION: 


In this text, we J ll have a a rigorous and an analytic treatment of the function of a 
vector state in space. We J ll build the concept using the properties of the following 
generalized topics: 


★ Linear Operators : That the V-l = i , -i* i = 1 and i 2 = -1 

■ That Polyphase AC is a Mathematical Model with NO INDEPENDENT COMPONENTS; only exists as SQUARES 
or PRODUCTS. 

■ Elementary Algebra of Complex Quantities: First., that the EFFECTIVE VALUE of a Sine Wave = 0. 

DC is represented by the following relation : Y = ■\f b 2 + g 2 and Z = ■\]r 2 + x 2 

Complex conjugate of AC is represented by the following relation : Y = g + ib and Z = r - ix 

Ohms Law is characterized by the relation E = IZ and I = EY where 

1 , , , 

Z = Impendance = — which is an absolute value derived from the equation Y = b + g and 

1 , , , 

Y = Admittance = — which is an absolute value derived from the equation Z = r + x 

11 11 

Proof of Z and Y : It is understood that the PRODUCT of Y 2 Z 2 = 1, so Z = — = - and Y = — = - 

Y V b 2 + g 2 z V r 2 + x 2 

E = Electromotive Force 
r = Resistance 
x = Reactance 
g = Conductance and 
b = Susceptance 


From the above: 

■ Y is defined as the resolving CURRENT into two COMPONENTS, and that it is in PHASE 
and quadrature with EMF, or the ENERGY CURRENT and the WATTLESS CURRENT, and 

■ Z is defined as the resolving EMF into two COMPONENTS, and that it is in PHASE and 
quadrature with CURRENT or the ENERGY EMF and the WATTLESS EMF 

In other words. Admittance is the reciprocal of IMPENDANCE, thatis. 



when we substitute for the COMPLEX CONJUGATE relation of AC 
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It is now necessary to descrie briefly; the notion of FORCE, also called the MAGNETOMOTIVE FORCE, 
and its application in INDUCTION MOTORS. For the relation T = mF0Sina, where 
T = Torque 
m = Constant 
F = Resultant MMF 
0 = resultant Magnetism and 

Sina = Angle of hysteric advance of PHASE or the Magnetic Cycle 

It is observed that the APPARENT EFFICIENCY of the motor is a FRACTION of the 

primary hysteric ENERGY and is < (1-S) Sina, where S represents SLIP or SYNCHRONISM. 

■ This is the base prerequisite knowledge, and sets a fine fundamental formalism for 
our purpose. 

I strongly recommend an adequate comprehension of C.P.Steinmetz's excellent book titled, 
f THE THEORY AND CALCULATION OF AC PHENOMENA'’ before you proceed, or you will likely and 
completely run out of runway with this text. 

■ An evaluation of the structure of o specif ied ALGORITHM, 

its relation to PATTERNS OF IRRATIONALS, INTEGER ARITHMETIC, 

and the METHODS of HIGH PRECISION MULTIPLICATION. This is Computer Science ! 

Let AC = 500 000 000 000, where AC = CB = Radius of A Circle 

B 

Then AD = AC - CD = AC - "sj AC 2 - 1 = 500000000000 - y/500 000 000 000 2 - 1 



Observe that CD < UNITY. Here we ' 11 use the FORMULA AD = AC - Sqrt [ (AC) 2 - l], 

to GENERATE a PATTERN OF IRRATIONALS called the CATALAN NUMBERS . The ALGORITHM for COMPUTING Catalan 
Numbers is described as a procedure of DIVISION with MULTIPLICATION SUBROUTINES. That is, 

AD = n T = en 31og2 + — =en LS!4!62S , where the EXPONENT of (n) is an IRRATIONAL NUMBER. Therefore, 

2 

we can conclude that for 4 multiplies of — digit numbers taking [e (n 2 ) ] time to compute, then 
4 nT . 

M) = nT = 8n + —-— representing a LINEAR TIME for addition. 

Its not so obvious at this point that an understanding 
of ALGORITHMS will be indispensable in the analysis of the FUNCTION 
X, = b m (mod X u ) , noted as the CONGRUENTIAL MULTIPLICATIVE GENERATOR in Wolfram 
Language. (*There cannot be a passive treatment of RANDOM NUMBER GENERATORS*) 

■ Assuming we're now adequately prepared and equipped with the formal background to 
analyze and comprehend the STRUCTURE of VECTOR STATE SPACE, we'11 learn that 
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REAL NUMBERS are represented as COMPLEX NUMBERS by showing in this instance if 

COEFFICIENT or COMPONENT in a COMPLEX VECTOR., then its COMPLEX CONJUGATE is 

1 - i ’l 2 2 1 

= — = — such that 
4 2 


1 

, V2~. 

1-i 

2 


and that. 


the INNER PRODUCT or MAGNITUDE of the COMPONENT 
1 1 

-= - and thus is related as the square of the the length of the VECTOR . 

V2" 2 

Representation of EIGEN VECTORS and EIGEN VALUES will be evaluated inclusive of the 

PAULI MATRICES ( ' J ® * ), and ( ® ) without which, 

the STRUCTURE and DEFINITION of STATE SPACE would break into oblivion. 

Finally, the OBJECTIVE will be to illustrate the TRUE and CLEAR meaning of the idea 
of ABSOLUTE SPACE, or the notion of CONSERVATION OF ENERGY . 


This is only possible MATHEMATICALLY, as we shall see, 

by applying the elegant FOURIER TRANSFORM PAIRS constructed from a fixed POISSON 1 s COMMUTATOR RELATION 
and HEISENBERGS 1 s representation of the matrix formulation of the way averages change with time. 


Ja e - ip (x-x ) =5 (x-x') 
fa ( e lx ( P ~ P ) = 6 (P - P ) where ; 


6 = Dirac delta function 
P = Momentum 
x = Position 
c = Speed of light 

e ip <x-x’) _ NORMALIZED or ORTHOGONAL BASIS VECTOR or HERMITIAN COMPLEX CONJUGATE 
It follows then, that a CLEAR neat equation emerges from applying 

ai It a (x) 

Schrodinger ' s WAVE EQUATION — = -c —— or the COMPLEX FUNCTION 

a t a x 

1 A (x) = (x - ct) merged with the HAMILTONIAN description of the MOTION OF A PARTICLE in one DIMENSION 
moving along a line and in terms of a VECTOR STATE SPACE , yielding the conjugate relation If 

Q'h . , a« 

-= x = c, then -= -P = 0. 

ap ax 


P > 


h 

< - 
2 


We ' 11 conclude by defining ABSOLUTE SPACE as < x 







THE DERIVATIVE OF AN IDENTITY_V2_Mathematica.nb • 8/10/18 


5 
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1. SPRING-PROBLEM APPROXIMATION 


relation between MASS and POTENTIAL 


Consider the HAMILTONIAN description of the particle 
moving under the influence of the Force exerted by the spring 

p 2 i 

Ti = — + — (Kx 2 ) where 
2m 2 ' 

P = Momentum 

m= Mass of the particle 

K= Spring Constant and 

x = Position of the particle 

It is understood also that 


T = — mV and Classically, 


p = mV 

For a SMOOTH POTENTIAL, 

we shall later in the text re - evaluate Momentum as described by the STATE VECTOR 

dV dP 

P = F = - — = — which is the IME DERIVATIVE of MOMENTUM and is equivalent to FORCE . 

dx at 

For now, let us specify the initial condition 
as [x (0), P (0) ] and there is no coupling to an external system. 

If the Energy E is a CONSERVED QUANTITY, 
the particle will follow a trajectory on a surface of CONSTANT ENERGY, 
or SMOOTH POTENTIAL , or INERTIA, given by the relation 
Kx 2 p 2 

-+ = 1 . We can then describe the EQUATIONS OF MOTION 

2 E 2mE 

as the relation between VELOCITY and MOMENTUM by differentiation where 

dx aw p 

— = -= - is the MOTION OF A WAVE PACKET or the relation between VELOCITY and MOMENTUM, 

at dP m 

dP /aw\ 

while — = - - = -Kx implies that the TIME DERIVATIVE of MOMENTUM is FORCE . 

at [dx) 

■ Suppose we would like to compute some PROPERTIES of the SYSTEM as it moves along a 


PATH 


[x (t), P (t) ] in PHASE SPACE using NUMERICAL 
INTEGRATION. (*We shall later establish the equations of 
motion in vector state space for an analytical treatment*) 
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We go about the solution of APPROXIMATING the CONTINUOUS PATH by a FIRST ORDER DIFFENTIAL taken as 

df 1 . . 

— = — (/(h + t) (*Refen to Euler's Method of Linear Approximation*) 


(h) being the BASIC TIME STEP. Here the solution of the equation 
of motion is obtained at times which are multiples of the basic unit (h) . 

It follows that a RECURSION FORMULAE for the position and momentum can be 
obtained by inserting the differential into the equation of motion as 

ax 1 p (t) 

— = — (x (h + t) - t x) = —- and 

dt h m 

dP 1 

— = — (p (h + t) - tP) = -Kx (t) or simply reshuffled as 

at h 

hP(t) 

x (h + t) = -+ (t) x and 

m 


P(h + t) = Pt - Khx(t) . 

So given an initial position x(0) and momentum P(0) consistent with a given Energy, 
the TRAJECTORY of the particle is SIMULATED at 
t = 2 h, 311, _, etc noting that the TOTAL ENERGY is set to UNITY. 

■ So, the trajectories in PHASE SPACE for the SPRING PROBLEM is thus GENERATED with a 
simple ALGORITHM where if the path is spiral indicates that the Energy is absorbed 
by the particle. A more reasonable trajectory is the result indicating an almost 

CONSERVED ENERGY. 

■ This method simulates a simple and elegant approach to the spring problem. We are 
therefore interested in developing ALGORITHMS OF HIGH ORDER in (h) to keep the 
ERROR as small as possible, and a TIME-STEP that the SYSTEM reaches an EQUILIBRIUM 
during a run. 


Suppose you were to Develop an algorithm and write a computer program for the pendulum 

d 2 <p i g \ 

-= - — Sin 0 where (1) is the length of the string and (g) is the gravitational constant. 

at 2 ' L I 

This is a SECOND ORDER differential equation, and the solution involves PARTIAL DERIVATIVES. We 
shall look into the method of INTEGRATION BY PARTS at a much later context. 


ORGANIZATION OF SIMULATION 

No matter how high the ORDER OF AN ALGORITHM, a system will eventually depart from the 
trace trajectory, as long as (h) is FINITE. 

A drift in the energy <5E is caused by the finite time step. This raises the 
question for the conservation properties of the algorithms during the simulation. To 
establish conservation is to CONSTRAIN the system the system artificially; 

Instead of using Forces to calculate the motion we use POTENTIALS. In addition, 
we may ask for the TIME REVERSAL PROPERTIES noting only the one step method is 
INVARIANT under time reversal if we require that the equations define a TRANSFORMATION. 

What may be the reason for Energy fluctuations? Why are the trajectories then NOT time- 
reversible? 
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To constrain the Energy, a FEEDBACK LOOP is established which introduces a TIME DELAY 
from the differential equations leading to fluctuations in the AVERAGE ENERGY(less signifi¬ 
cant though). 

It is desirable to have an algorithm which changes a value from step to step. Initially, 
the values should be large to compensate negative results, and gradually decreases as 
the system approaches equilibration. What would be a computational scheme for the TRUN¬ 
CATED OCTAHEDRON boundary condition? 


DEMONS 


Let us now develop an algorithm to specify the TRANSITION PROBABILITY from one state (x) of the 
system to a state (x') per unit time. The goal is to ensure positivity of the domains's Energy. 
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First, we note the transition probability as [W (x, x ') ] . To ensure that the states are ultimately 
DISTRIBUTED according to [P (x) ] equilibrium states, restrictions must be placed on [W (x, x')]. 

For all complementary pairs (s, s) of sets of PHASE POINTS, 
there exists (x 3 S) and (x3s] such that [W (x, x ') t 0] . This means, for all (x, x '), [w (x, x ') > 0 ]; 
for all (x) the k/(x, x ') =1; and 

for all (x) the w(x, x ') P(x) = P(x) . 

The DISTRIBUTION of STATES must corespond to the EQUILIBRIUM DISTRIBUTION [Z’ 1 f [H, {x} ] ] 
because [P (x) ] = [z _1 f[H, {x> ] ] where (Z) is the proportionality constant and does 
not enter transition probabilities noting also its not accessible in a simulation. 

This implies that the free ENERGY [F = -K B Tin (Z) ] cannot be computed directly. We must then 
compute properties of the system using the proportionality constant (Z) since the KINETIC 
ENERGY term from the HAMILTONIAN is dropped and thus we cannot use the equations of motion . 

Now, suppose that a state (x) is generated such that [H (x) ] E representing a RANDOM 
WALK ON A CONSTANT SURFACE IN PHASE SPACE. A sampling algorithm on the surface then 
has to produce further states and allow for (3 ) variations or equilibrium in the 
region [E - 3 < H (x) < E + 3] away from the surface. We do this by introducing a DEMON, 
known as a EXTRA DEGREE OF FREEDOM with Energy (E D ) into the proportionality constant such that 
Z = <5 [H ( x ) + E d - E] where the distribution is 

expressed by a DELTA FUNCTION. So the AVERAGE VALUE of Z = J~<5[H(x) - f] dx. 

You must know that With any OBSERVABLE (A) is associated a FUNCTION [A (x) ] which 
is the MEAN AVERAGE.(* Refer to Probability Density Function *) . Thus, for 

A = - f A (x) 6 [H (x) - E] dx 

2 Ja 

implies sampling the available PHASE SPACE 
of the system [a] and carrying out a SUMMATION of the function (Z) . 

The DEMON plays a role similar to the Kinetic Energy by producing changes in the 
configuration by traveling around the system and transferring Energy. It 
thus creates a Random Walk of the system on the surface. However, 
we MUST restrict the Demon ' s by limiting its Energy to positive values otherwise 
it will absorb all the energy (compare with step 7 in the algorithm below) . 

Algorithmically, an outlined and summary of the procedure takes the form: 

1. Construct a state such that [H (x) = E] 

2. Set the Demon Energy (E D ) for example [ (E D ) = 0] 

3. Choose a part of the system 

4. Change the local state of the system so that (x) truncates to (x ') 

5. Calculate the Energy change produced for instance [AH = H (x ') - H (x) ] 

6. If the Energy is LOWERED, accept the change, 

set E d <— (E d - A h ) and count (x ') as a new configuration. Return to step 3. 

7. Otherwise, accept the change ONLY IF the Demon carries enough Energy, 

that is (AH - E d ) > 0. In this case, E D <— ( E D - A H ) and count (x ') as a new configuration. 

8. Return to step 3 

The algorithm guarantees with step 6 and step 7 that the system relaxes to thermal equilib¬ 
rium. In addition, step 7 also ensures the positivity of the DOMAIN'S energy. 

Conceptually, we may view the Demon as a thermometer. Indeed, the Demon can take up or 
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lose Energy as it is successively brought into contact with the parts of the system. 

Initially, the Demon has an arbitrary distribution. The system thus acts as a reservoir 
and renormalizes the Demon. Ultimately, the energies become BOLTZMANN DISTRIBUTED, allow¬ 
ing the calculation of the temperatures as 

P (E d ) oc exp (- 

V T Kg 1 

Below is a schematic representation of a random walk on 
a constant Energy surface is Phase Space using the Wolfram Language 
(*SeedRandom[1234] ;data2d=RandomFunction [RandomWalkProcess [0.5], {0,10*4},2 ]; 

Graphics[Line[Transpose@data2d["States"]],AspectRatio->Automatic]*) 



■ We can see more clearly the meaning of the choice of transition of probabilities 
where the system is driven towards the minimum Energy. Thus, configurations which 
raise the energy are only accepted with a Boltzmann Probability. 


LOGICAL OPERATORS 

The self-avoiding random walk is where the walker does not cross its own path. At each 
step, the walker checks if the neighborhood sites have been visited before. Of course, 
the walker is not allowed to retrace its steps. Quite often, the walker encounters a 
situation where all the sites in the immediate neighborhood have been visited before; The 
walk then terminates. 

To be able to write a program which shows a screen how a self-avoiding random walk pro¬ 
ceeds across a two-dimensional grid, one has to assume a RECURSIVELY SUMMING PSEUDORANDOM 

NUMBERS. 

■ Before we proceed further in the analysis of RANDOM NUMBER GENERATORS, get a clear 
understanding of the methods of LINEAR CONGUENTIAL GENERATORS and FEEDBACK SHIFT 
REGISTERS. 

■ Run the programs to get acquainted with the workings of the algorithms and to 
develop the requirements of a simulation. 

■ Find a thorough treatment on how to generate RANDOM NUMBERS in Donald Knutlr’s f The 
Art Of Computer Programming Volume 2 (Chapter on Semi-numerical Algorithms). Hands- 
on experience with the C-Programming Language expected. 
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■ Note that such an algorithm is an ITERATION of a fixed LINEAR TRANSFORMATION on a 
FIXED VECTOR SPACE. 

By a general THEORY of LINEAR RECURRENCE, each entry in the ARRAY(VECTOR SPACE) is a 
LINEAR RECURRENCE SEQUENCE corresponding to the characteristic POLYNOMIAL. 


For [(n*w) - r] array with a characteristic polynomial [0B (t) ] of the TRANSFORMATION (B) , 
the sequence attains the maximal PERIOD [2 P - l] = 2 <nw_r) - 1, 
if and only if [0B (t) ] is PRIMITVE (has ROOTS) 

The test of PRIMALITY of an INTEGER is achieved by eliminating (r) BITS from 
the (n*w) array so that the DIMENSION of the STATE SPACE is (nw-r) array. That is, 

PRIMES of the form (2 P - l) upto (P = 1398 269) need no FACTORIZATION .We ' 11 later 
understand the implication of COMPONENT FACTORIZATION and the concept of CORRELATION. 

In determining the next STATE, each BIT must be fully REFLECTED thus results the recurrence, 
and by MULTIPLYING with a MATRIX (T) , called TEMPERING, a good K - Distribution is realized. 
(*To fully appreciate the concept, get acquainted with the 
Multiplicative Recursive Matrix Method and Inverse Decimation Method*) 

Its interesting to note that in the case of (P) being very large ( > 10000), 
the direct computation may need several years to catch a primitive 
polynomial. A concise illustration is in the application of TIME COMPLEXITY. 
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2. AN INTRODUCTION TO CRYPTOGRAPHY 


For IT Security, its relevant to understand the classification, basic setups, SUBSTITU¬ 
TION CIPHERS and attacks on crypto-systems. Some of the modern applications of Cryptogra¬ 
phy include: 

■ Cell Phones (GSM Voice Encryption) 

■ Bank Cards and Online Banking 

■ Electronic Passports 

■ Secure Shell, E-mail Encryption and Plug-ins for Web Browsers 

These applications entail a two diverse unified subjects of CRYPTOLOGY. That is: 

■ Cryptography; Involves Symmetric, Asymmetric Algorithms(PKI) , and Protocols 

■ Crypto-analysis 

The setup for symmetric cryptography raises the problems of communication over an inse¬ 
cure channel. Here the channel examples are the INTERNET, and WI-FI. 

Note that Public Encryption Functions are the only ways so far, of keeping these channels 
SECURE(defined as CRYPTOANALYSIS). 

In practice, one should never use an untested CRYPTO-ALGORITHM. In principle, a CRYPTO¬ 
SYSTEM should be secure even if the attacker knows all the details about the system, with 
the exception of the SECRET KEY. Observe that the principle is counterintuitive. 

The SUBSTITUTION CIPHER is a historical cipher that operated on letters of the alphabet. 
For example, if 

A represents m 

B represents d 

C represents w 

E represents q and so forth....We obtain a fixed look-up table. For instance, e[CBBC] 
represents (wddw). Is this cipher secure? How can it be broken? 
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1 . 


2 88 divided into (2 56 * 2 32 ) resulting in = 

10 9 days TIME COMPLEXITY. This implies the keyspace it too large. 

2. Letter Frequency Analysis where the letters have very 
different frequencies in the English Language. Moreover, 
the frequency of PLAIN TEXT letters is preserved in the CIPHER TEXT. 

For instance, the letter (e) is the most common letter in English. Almost 13 % 
of all letters in a typical English text are (e) . The next most common letter is (t), 
with about 9 %. This method of analysis works because identical plain texts map to identical 
cipher text symbols thus Its weakness.So if an attack is succesful, then the system collapses. 


There are often many possible attack approaches used to describe cryptosystem attacks. 
These are: 

1. Classical Crypto-analysis; Includes the briefly described Brute-force and Analytical 
attacks(letter frequency analysis) 

2. Social Engineering 

3. Implementation attacks, where the HARDWARE is not protected. This is an interesting 
“side Channel analysis” of the key SPACE, where one records the power TRACE (walk) on an 
oscilloscope once the algorithm on the microprocessor has been triggered. The timing 
behavior (step) depends on the key BITS (array). We then use the timing code to com¬ 
pute the key. 


DEFINITION : 

A PSEUDORANDOM SEQUENCE (Xi) of (w bit) integers of Period (P), 
satisfying the following conditions is said to be K - Distributed 
(*Refer to our discussion of our discussion on TEMPERING*)to (v_bit) accuracy : 

Let [Trunc v (x) ] denote the number formed by the leading 
(v_bits) of (Xi) and consider (P) of the (Kv_bit) vectors such that 

[{Trunc v (Xi), Trunc v (x 1+1 ), ., Trunc v (x 1+( k.i))}] for [0«i<P] 

where [ (i= 1, 1, ... .P - l) ] represents coordinates and 
(K) represents the DIMENSIONAL UNIT 

The SEQUENCE is said to be K - Distributed to (v_bit) accuracy if each DIMENSION contains 
the same number of points (except for the origin) . Consequently, the higher (Kv_bit) 
for each (v_bit) assures a higher dimensional EQUIDISTRIBUTION with (v_bit) PRECISION, 
based on a LINEAR RECURSION over the two ELEMENT FIELD (F 2 ) , 
known as the (F 2 ) GENERATORS.(*This should be apparent to you, 

considering that a linear sequence has a characteristc polynomial*) 
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The cryptographic interpretation of the K-Distribution will assure that such a sequence 
is K-Distributed to (vbit) accuracy, and that all the SEED(INITIAL) BITS are RANDOMLY 
GIVEN. Thus, the knowledge of the most significant (v_bits) of the first (w) words does 
not allow the user to make any statement about the most significant (v_bits) of the next 
word, if (w < K).This is because every BIT PATTERN occurs equally likely in the (v_bits) 
of the next word, of course, by definition of K-Distribution. 

Thus, if the simulated system is sensitive only to the history of (K) or less previously 
generated words with (v_bit) accuracy, then it is THEORICALLY SAFE. (Refer to a detailed 
justification of the recurrence in the paper “MERSENNE TWISTER: a 623-Dimensionally 
Equidistributed PseudoRandom Number Generator” by Makoto Matsumoto and Takuji Nishimura) 


Another and important criterion criterion on the RANDOMNESS of the (F 2 ) GENERATORS is 
in the number of terms in the characteristic polynomial of the STATE TRANSITION FUNCTION 
(♦characteristic polynomial with many terms*) 


We J ll later compare the speed (CPU TIME) and WORKING AREA of pseudorandom number genera¬ 
tors; Although high speed is probably due to the hardware architecture like CACHE MEMORY 
and PIPELINE PROCESSING. 


■ Note that the calculation of the recurrence is realized with BITSHIFT. 
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3. PATTERNS OF IRRATIONALS, AND THE COMPLEXITY OF AN ALGORITHM 
Its understood that the ratio of the diagonal on the square is INCOMMENSURABLE. 


An output of the (V2 j in Wolfram Language, N[Sqrt[2], 18] outputs 1.41421356237309505 
We now can define a Catalan Number as the 

Set (P) of balanced parenthesis strings where [A3P] where A is the empty string. Thus, 

If [(a, /3) 3P], then [(a) /3 3 P] 

We can obtain a balanced parenthesis string via the second rule, from a unique (a, /3) pair. Therefore, 
an ENUMERATION of the CARDINALITY OF THE SET can be constructed as follows : 

Let (C n ) be the number of balanced paranthesis strings with exactly (n) pairs of paranthesis, 
where C 0 = 1 and is the empty string. Then C 2 = 1 => ( ), 

C 2 = C 0 Ci + C, C 0 where we choose the a, then /3 ; Thus C 2 = 2 ; and 

C n = ^ C k C„_ fc n > 0 representing the GENERATOR. 

/f=0 


Consider the Newton ' s Method of approximating SQUARE ROOTS 
where {y = f(x) } . Then the derivative j y' = = SLOPE OF THE TANGENT LINE 


We find the roots of/(0) = 1 through succesive approximation. For example, 
we take the equation [y=/(x)] and note the TANGENT to this curve. Thus, 


the equation of the TANGENT is the DERIVATIVE and so, 
x 0 = 1.000000000; where we start with one digit PRECISION 
Xj = 1.500 000 000 
x 2 = 1.416 666 666 «-» 2 

x 3 = 1.41215686 <-> 4 

x 4 = 1.414213562 <-> 8 


The important observations in this recurrence is that. 

Integers with certain digits have a QUADRATIC CONVERGENCE where the number of digits double 
with every ITERATION resulting in a LOGARITHMIC number of iterations. For instance, 

note that (2, 4, 8, 16, _, ...) is a GEOMETRIC SERIES. 

If we take (a) to be an integer or fraction or a small number, and (Xi) to be a number with large digits, 
we need COMPUTERS to divide (100 _ digit numbers) by (50 _ digit numbers) .So, 

for [ (—) * (—) ] is referrred to as a HIGH PRECISION MULTIPLICATION 


because it is a division algorithm with multiplication subroutines. 

If we assume (V2~j to be a double digit precision as from above, 

then we want an integer that is j Tp2 10 d j = |a/2 x 10 2d j (*The symbol '==' 


means Tilde_Full_Equal*) 


This is the definition of the algorithm for applying 
Newton ' s Method to form the desired high precision multiplication . 

Note that (d) is the number of DIGIT PRECISION. For instance we may request for (d = 10 6 ) . 


Let two (n_digit) with Radix (r = 2, 10) where the RANGE [ 0< (x, y) < r 11 ] is precise. 

■ We J ll use the strategy of DIVIDE and CONQUER where we break the digits to 64_BIT, 
and let the machine do the rest of the computation. 
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Let { (x = Xi r n/2 + x 0 ) }; and 

{ (y = Yi r n/2 + ye) } Where (x 2 ) is the high half and (yi) 
is the lower half. (* Obserser the characteristic polynomial*)Then, 

{ (0 « (x 0 , x a ) < r n/2 ) } and {(0< (y 0 , yi) « r n/2 ) } thus a straighforward 
decomposition of the multiplication operation; 
let { (z 0 = (x 0 y 0 ) }, 

{ (Zi = (x 0 y 2 + x 1 y 0 ) } and thus, 

{(z 2 = x 2 y 2 )} (*Refer to the enumeration of cardinality of the set, 
at the beginning of the chapter*) Then, (z = (xy) = x 2 y 2 r n + (x 0 y! + X!y 0 ) r n/2 + x 0 y 0 )j 

(*We've applied and substituted the characteristic polynomial 
function to obtain the formula of the algorithm*) 

Observe now that the (z) has 4 multiplies of (n / 2) digit numbers at the second term [ (x 0 y ± + Xx y 0 ) r n/2 ]; 
which is equivalent to 8 (n 2 ) TIME to compute (* The log 2 4 = 2 *) .So 
we can assuredly remark that z) is a e (n 2 ) ALGORITHM. 

Then in GENERAL, T (n) = 4 T (n / 2) + e (n) represents a LINEAR TIME for ADDITION, 
calle d the COMPLEXITY OF THE ALGORITHM. 

We can reduce the complexity further using an algorithm with 3 multications instead of as of 0 (n 2 ) . So, 
If T (n) 3T (n/2) +6 (n) =e (n log23 ) =e (n 1 - 5848625 ) t 
meaning the complexity of the algorithm is an IRRATIONAL NUMBER. 

At the beginning of our Introduction, we supposed a really big circle with center (C) and 
a Diameter of 1 Trillion. We then dropped a perpendicular one UNIT high so that 


CB = CA = r thus AD = AC - CD = 500 000 000 000 - 500 000 000 000 2 - 1 



CD is not a perfect square and so is an irrational quantity. So by sutracting 1 
from it we obtain an irrational quantity that goes to infinity . The algorithm for 


AD = AC - CD computes the CATALAN NUMBERS. 
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4. COMPUTATION COMPLEXITY 

The MERSENNE TWISTER operates with (P = 19 937) which is a PRIMITIVE twisted 
generalized FEEDBACK SHIFT REGISTER SEQUENCE (*Review self-avoiding random walk*) 

The algorithm does not involve any arithmetic since all operations are SHIFTS, 
AND's,OFPs,XOFPs. 

A SEED(single integer) initiates the whole process generated upon start up with an 
unsigned 32_bit random integers. 

There is TEMPERING with logical operations to improve RANDOMNESS. So, when the pointer 
reaches the end of the cache where the last element of the STATE is a pointer in cache, 
its refilled with another 623 ELEMENTS. 

Our purpose is to analyze the group theoretic properties of PERMUTATION and TEMPERING 
operations. 

First, we note that the PERIOD of the Mersenne Prime is (2 19937 -1). Its been demonstrated 
that the design generates results that satisfy an EQUIDISTRIBUTION property in 623_Dimen- 
sional cube. (Review the C_code for MT19937) 

The code works with both 32_bit and 64_bit machines. 

The function [genrand () ] returns a UNIFORMLY DISTRIBUTED REAL PSEUDORANDOM 
NUMBER of TYPE DOUBLE with 32 bit PRECISION in the closed interval [0, 1] . 

The function [sgenrand ()] sets the INITIAL VALUE to the ARRAY {mt [N] } .Before 
using [genrand () ], [sgenrand () ] must be called with a NON_ZERO unsigned long integer as SEED. 

There are '•’Magic Numbers” in a macro table for replacing FUNCTION CALLS to enhance RUN¬ 
TIME. The GENERATOR has MULTIPLIERS that attribute a LINEAR CONGRUENTIAL SEQUENCE with 
extremely LOW SERIAL CORRELATION over the entire PERIOD, where the ACCURACY of the RANDOM 
SEQUENCE should be the SAME IN ALL DIMENSIONS. 
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RANDOM NUMBERS are truncated say to (lgV t _bits) where (V t ) decreases with increasing 

(t) and the (t_DIMENSION) is the the source of the random number |-\J x 2 2 + .... +x 2 t J 

such that (x a + .... +a t-1 x t ) s 0 (modulo m) . (*A low dimension 
being the most damaging to randomness*) 

The most generalized REAL VALUED DISTRIBUTION can be expressed in terms of its DISTRIBUTION 
FUNCTION/(x) , which specifies the PROBABILTY that o random quantity (X) will not exceed ( x ); 
Thatis/(x) = P(r) [X « xj always increases monotonically from 0 to 1. This means/(X!) < 

/(x 2 ) if x a < x 2 and where/ (-oo) =0and/(+oo) =1. 

If/(x) is continuous and strictly increasing so that/(x a ) «/(x 2 ) whenx!<x 2 , 

it takes on all values between 0 and 1 and there is an INVERSE FUNCTION/ -1 (y) such that, for 0 < y < 1, 
y = / (x) if and only if x = / -1 (y) . 

We can thus compute a random quantity (X) with distribution 
/(x) by setting (X) = / -1 (u) where (u) is UNIFORM with the PROBABILITY that (u) «F(x) . 

It follows then that NUMERICAL ANALYSIS is considered as good 
method for evaluating/ -1 (u) to a desired accuracy. 

Understand that any algorithm that employs a random number as input will give a random 
quantity with some distribution as output. The task is to find general methods for con¬ 
structing the algorithm, given the distribution function of the output. 

From a Theoretical point of view, we consider optimal ways to generate random variables 
with a given distribution in the sense that the method produces the desired result from 
MINIMUM possible number of random BITS. For instance, in the HARD_CORE_BIT method, we 
consider 

a Boolean function/such that/(x) and g (x) are easily computed, 

although/[g -1 (x) ] isNOT.(*For large modulus generators, the primitivity check 

is a hard number theoritic task to ensure a reliable source of random numbers*) 

This is the reason for having a linear congruential sequence with low serial correlation 
over the entire period. 

In PROGRAMMING, the absolute difference between two expressions could be evaluated at 
randomly generated points to test for INEQUALITY OF THE EXPRESSION. 

For instance, there is no evidence that the yj x 2 and the | x | are different for REAL VALUES; 
while there is evidence that the x 2 and the | x | differ for at least some COMPLEX VALUES. 

The point is, there are features in Programming Languages that seem a bad fit for the 
language. These includes the brevity of the language caused in part by the excessive 
reuse of symbols and problems with OPERATOR PRECEDENCE. Worse than that, many symbols and 
keywords are “overloaded”; given different meanings when used in different contexts. 

These make the SCOPE of the RULES not intuitively clear. 
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For example, the C symbol " () " has multiple different meanings 

1. Enclose formal parameters in a function definition 

2 . make a function call 

3. Provide expression precedence 

4. Convert (cast) a value to a different type 

5. Define a macro with arguments 

6 . Make a macro call with arguments 

7. Enclose the operand of "size_of" operator when it is type name 
So, What does, "apple=size_of(int)*p" 

mean ?Is it the size of an "int" multiplied by "p" ?0r the size of whatever "p" points but cast 
to an int ?0r something even weirder ?Can one write a program to probe questions like this ? 


The more work you make one symbol or expression do, the harder it is for a compiler to 
detect anomalies in your use of it. This is an example of the C_Programming Language a 
little further out on the ragged edge of ambiguity! Further, it is not surprising that 
any change in PRECEDENCE will impose an intolerable burden on the existing source base 
and importantly the order of operand evaluation. 

The solution is to always put parenthesis around an expression that mixes Booleans, arith¬ 
metic,or BIT_TWIDDLING with anything else. Besides, why is the order in which these group¬ 
ings are evaluated undefined? 

We J ll now illustrate excellent and explicit application of parenthesis around expres¬ 
sions . 

Below is a "random walk” generated in MATHEMATICA by recursively summing pseudorandom 
numbers, starting at 0; 

ListLinePlot[loin[{0. }, Accumulate[RandomReal[{-1, 1), {100}]]]] 

Out = 



Test for the equivalence of symbolic expressions by substitution of random numbers; 
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Max [Abs [Sqrt [x A 2] - Abs [x] / . x -» RandomReal [ {-10, 10}, 10 000] ] ] 

Out = 0. (* This provides NO evidence that 
the Vx 2 and the|x| are different for REAL VALUES *) 

Max [Abs [Sqrt [x A 2] - Abs [x] / . x -> RandomComplex [ {-10 + 10 I* 10 + 10 1}, 10 000] ] ] 
0ut=14.142(* This provides evidence that the 
Vx 2 and the |x| differ for at least some COMPLEX VALUES *) 


DEFINING RANDOM NUMBER GENERATORS IN WOLFRAM LANGUAGE :Writing the Code 

MULTITIPLICATIVE CONGRUENTIAL GENERATOR 

(*A multiplicative congruential generator follows the RECURRENCE RELATION [Xi=b (x^) mod m]*) 
(*The generator defined below, will alow ONLY for generation of REAL NUMBERS*) 
(♦Initialization of the generator will extract the values of the multiplier and 
modulus.Initialization will fail if either of these values is not a positive integer*) 
(♦Calls from the kernel to Random InitializeGenerator are effectively wrapped in 
Catch.Throw is used in the initialization code to easily exit in case of problems*) 

(♦"True" Establishes that MultiplicativeCongruential generates REALS*) 

(♦The following initializes the generator*) 


Options [MultiplicativeCongruential] = ["Multiplier" -» 123456789, "Modulus" -» 2 A 35 - l]; 

MultiplicativeCongruential /: Random'InitializeGenerator[MultiplicativeCongruential, opts_] 

Module[[mult, mod, flops = Flatten[[opts. Options[MultiplicativeCongruential]}]}, 
mult = "Multiplier" /. flops; 

If[i [lntegerQ[mult] && Positive[mult]), 

Throw[$Failed] 

b 


mod = "Modulus" /. flops; 

If[l [lntegerQ[mod] && Positive[mult]), 

Throw[$Failed] 

b 

MultiplicativeCongruential[mult, mod, 1]]; 

MultiplicativeCongruential[_] ["GeneratesRealsQ"] := True; 

MultiplicativeCongruential[mult_, mod_,_] ["SeedGenerator"[seed_] ] : = 

MultiplicativeCongruential [mult, mod. Mod [ (mult* seed), mod] 
MultiplicativeCongruential[mult_, mod_, s_]["GenerateReals"[n_, (a_, b_), prec_] ] : = 
Module[{x = s), 

[a + (b-a) Table[x = mult*x; 

Mod[x, mod], (n}]/mod, MultiplicativeCongruential [mult, mod, x] ] 

1 

BlockRandom[ 

SeedRandom[Method -> MultiplicativeCongruential]; 

RandomReal[{15, 55), 5] ] 

(♦Generates 5 reals using the MultiplicativeCongruential generator*) 

Out : {35.19237603468915, 42.523274676278334, 49.298336948102175, 

42.65275420468687, 16.116889654691242} 
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BlockRandom[ 

SeedRandom[Method -> MultiplicativeCongnuential ]; 

RandomReal[{15, 55}, 5]]j// Timing 
Out : {0., Null} 

(* Add "//Timing" to see that the 

built in methods take almost no TIME TO COMPUTE the values *) 

BlockRandom[ 

SeedRandom[Method -» MultiplicativeCongnuential ]; 

RandomInteger[{5, 50}]] 

Out : 

... Randomlntegen : "The current random generator does not support generation of random 
integers in the given range." 

Randomlnteger[{5, 50}] 

(* It means the generator is not defined for integers *) 

A Random Distribution Vector can also be used to define generators for multidimensional 
distributions. 
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Lets suppose a random point from a UNIFORM DISTRIBUTION on the UNIT DISK where, 
the set of REAL points {x, y} with x 2 + y 2 < 1 

is desired. Such a RANDOM point can be constructed as follows : 

(* Generate a RANDOM ANGLE e uniformly distributed on [0,27r] *) 

(* Generate a RANDOM VECTOR (u) uniformly distributed on [0,1] *) 

(* Return jsine Vu”,Cose VtTJ. The returned ORDERED PAIRS can be MULTIPLIED 

by (r) to generate points uniformly distributed on a disk of radius (r) *) 

Below is the Wolfram Language Code : 

(* "RandomReal" defines a generator for a uniform disk of radius (r) *) 

(* "ListPlot" VISUALIZES the DISTRIBUTION of (l0 4 ) generated points on this disk *) 

UniformDisk/: 

RandomDistributionVector [UniformDisk [r_ ? Positive], n_Integer, prec_ ? Positive] : = 
r* Sqrt [RandomReal [l, n, WorkingPrecision -» prec] ] * 

Transpose [ [Cos [tt], Sin [tt] ]&[RandomReal [2 Pi, n, WorkingPrecision -> prec] ] ] 

RandomReal[UniformDisk[2]] 

Out : {-0.5549869088775746, 0.048327490357459944} 

ListPlot[RandomReal[UniformDisk[2], 1000], AspectRatio->l] 

Out : 
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5. VECTOR STATE SPACE 

DINSTICTION OF THE STATE OF A SYSTEM 

The AVERAGE VALUE of a VECTOR ct with an ANGLE of Cose for a q bit is + 1 or -1. 

Note that the STATE is NOT a SET. This makes a VECTOR STATE a collection of MATHEMATICAL OB1ECTS, 
that is COMPLETELY ABSTRACT VECTORS. But first, keep in mind the notion of a COMPLEX CONJUGATE. 

A COMPLEX NUMBER is a COMPLEX VECTOR SPACE. 

FUNCTIONS can also form a VECTOR SPACE. 

We shall RESERVE the terminology of a square matrix for it doesnt apply in this case. 

A COMPLEX CONDUGATE VECTOR SPACE, is a SPACE of COMPLEX CONDUGATE VECTORS. 

Take the symbol "<>" to represent ANGLULAR BRA_KETS. 

The Rules of Arithmetic Operations shall apply in all cases henceforth . 

Let Z | a > «-► < a | Z* 

If | a > implies j j, then < a | implies (a!*, a 2 *) 

(* Laid in a ROW to reprsent a complex conjugate vector space *) 

The INNER PRODUCT or DOT PRODUCT | a > <b| = <b|a>and note that <b|a> + <a|b> 
although they relate to the complex conjugates 
So, < b | a > = < b | a >* ; 

(/3 i*/32*) ( q 1 ) = Ai* “i + /3* 2 <*2 (* This is the sum of the product the COMPONENTS*) and 

(a* i«* 2 ) (^) = a*! jBi + a* 2 0 2 (* The complex conjugate of the above *) 

This is an illustration of a complex vector space . 

The INNER PRODUCT of a VECTOR <a|a> = <a |a> = REAL and + ve. Thus, 

(a*i a* 2 ) | j = a*! a 2 + a* 2 a 2 = REAL and + ve. 

This is equivalent to the SQUARE of the VECTOR. 

It is obvious and not obsecure to say that complex vector spaces have REAL and + ve LENGTHS. 

If the INNER PRODUCT < b | a >= 0, It is DEFINED as being iTHOGOf or PERPENDICULAR . 

The DIMENSION OF A VECTOR SPACE is the maximum number of ORTHOGONAL VECTORS you can find 
in the vector space. Also referred to as MUTUAL ORTHOGONAL VECTORS. The reverse of the 
statement is TRUE. 
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For instance, ( * ) ( ® j are ORTHOGONAL VECTORS because 1 (0) + 0 (1) = 0. 

This is a 2 DIMENSIONAL VECTOR SPACE and there cannot be a THIRD DIMENSION because of MUTUAL COMPONENTS . 
REAL NUMBERS are a special case of COMPLEX NUMBERS and that, (®) is ORTHOGONAL TO EVERYTHING. 


The OUTER PRODUCT is represented as | o> <b | (* Observe the NOTATION *) 

It is important to note that COMPLEX FUCTIONS form COMPLEX VECTOR 

SPACES where the VARIABLES are COMPLETELY UNCORRELATED. (* RANDOMNESS *) 

The AVERAGE VALUE of < o m > = n * m 

which implies the Cosine of the angle between (m) and (n) or the COMPONENT of (n) along (m) . 

The RULES OF CALCULUS OF PROPOSITION LOGIC (AND) & (INCLUSIVE OR) apply. Let, 

A : ct z = +1 
B : a x = +1 
while Z = +1 , 

A or B is TRUE, B orA is TRUE 
Then, 

< A | B > = COMPLEX VECTOR = COMPLEX CONJUGATE OF < B | A > . We also know that | A > + | B > is a KET VECTOR; 
So Z | A > is KET VECTOR multiplied by a COMPLEX NUMBER. 

Recall that < A | A > = REAL , 

which is analogous to the SQUARE OF THE LENGTH OF A VECTOR 
(* That is, the INNER PRODUCT OF VECTOR and ITSELF *) 

< A | B > =0 is an ORTHOGONAL (MUTUAL) VECTOR SPACE which is now an UNAMBIGOUSLY DISTINGUISHABLE 
(* That is, the INNER PRODUCT OF TWO DIFFERENT VECTORS *) 

These are the RULES for distinguishing a VECTOR SPACE. 

Lets now examine the PROBABILITIES OF VECTOR SPACES. Disregard whether the coordinate 
system is right or left handed. 


Note that {-i*i = lj and the V-l = i; 

So {e ie } will be considered as a PHASE FACTOR or UNIT LENGHT 
Consider the following VECTOR STATES in a coordinate system where 

I R > 

I L > 

| I > and 

|0> represent RIGHT, LEFT, IN, and OUT while 
| u > and 

| d > represent UP and DOWN 

For the KET VECTORS is a state of | u > and | d > , 

< u | d > = 0 

Then the PROBABILITY (sum of square of the length) for 

I R> = —==: ( | u> + | d>) 

V2 

I L > = ( | u > - | d >) 

V2 

(* The probabilities are the MAGNITUDES 
of the VECTORS or the square of the length of a vector *) 
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Thus < R | L > = 0 which is analogous to complex conjugation 
1 


I 1 > = 
I 0 > = 


VT 

i 

VT 


u > + 


u > 


VT 

i 

VT 


d > and 

d > 


If 


We then consider COMPONENTS for multiplying vector spaces 
0 ' 


■>-C) 


and d > 


-C)< 


Take |U> = |L> + |R> and |d> = | L > - | R >; Then from above, 

I R > = 


1 

V2 

1 

> 1 L> = 

1 

1 

, 11> = 

1 

V2 

i 

, and | 0 > = 

1 

i 

k V2~ - 


vr 


. *\ pr , 


k VT j 


You can easily prove that | I> and | 0> 

are ORTHOGONAL since the rule for calculating inner products requires that <1 | 0> =0 
and also show that the row vectors have complex conjugates , 

1 

meaning the SQUARE OF THE MAGNITUDE OF THE COEFFICIENTS MUST = 

Thus, in the | I> STATE, the probability of | u > or | d> is most strikingly —. Here is the PROOF; 
Lets us illustrate | I > in terms of | L > and | R > . For 


I I> 


VT 


( | L > + | R >) = — ( | L > + ] R >) = - L>+- R>+- 


- - I 

” 2 I 

1-1 


L > + 


1 
2 

1 + i 


i 

|R> 

1 L> 

+ VT 

VT 

VT 

i 

i i 


R > + — 
2 

R > - - 
1 2 

l L > 


R > 


Its important to VERIFY that the 

1 - i 1 

complex number-multiplied by its complex conjugate MUST = —.Thus, 

f 1 - 1 ) ( 1 + i ) _ (l-i + i + l) 2 1 

4 4 4 2 ‘ 

The statement above applies to ALL the PROBABILITIES of STATES of other combinations. The states 
are SYMMETRICALLY RELATED to each other. Its an easy concept but very unfamiliar indeed. 


2D SPACE can now, clearly be represented as 2D VECTOR SPACE. A thorough understanding of 
TENSOR CALCULUS is required to proceed. Interestingly, it takes two REAL NUMBERS to pick 
a direction. 


















26 


THE DERIVATIVE OF AN IDENTITY_V2_Mathematica.nb ■ 8/10/18 


NOTE that when | d > and | L> are OVERLAPPING, 

the PROTECTION PROPOSITIONS are NOT COMMUNITATIVE. (* We’ll illustrate this further in the text *) 
Moving forward. 

If | A > = | i > then, 

1 

< J | A > = < j | i > and can be represented as <5 3i . This is the KRONECKER DELTA symbol. 

i 

So, the KET |A>=^|i> <i|A> where 

i 

| A > is written as the sum of the inner product of the BASIS or UNIT VECTOR and thus, 

the BRA < A \ = ^ < A | i > < i | . This is the EXPRESSION of the VECTOR SPACE in terms of its COMPONENTS . 
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6. LINEAR OPERATORS 
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An operator can be defined as a procedure (M) acting on | A > to give | B > . That is M | A > = | B > . 
SoM [Z | A >] =z M | A> (* The relation is equivant *) 

We ' re now intersted in the components of the output vector | B > . 

We do this by projecting the inner product of both vectors using (i); 

Thus for the case of of M | A > = | B > we obtain < i | M | A > = <i|B> = I3 ± . 

The ^ < i | M | j > < j | A > = < i | M | j > a-j = /3 { . So ^ Ai ± j cxj = /3 t . This can be illustrated 
i 

by multiplying a rule with a MATRIX by a column vector. For instance, 


Mu 

AJ 21 


We can the define < A | M | B > by noting [ < A | M] | B > = < A | [M | B >] . So now, 

the ^ < A | i > < t | M | j > <j|s> = . a*i Mij 

ij 

If M | A> = | B >, 

we note that <A | M = <B | is COMPLEX CONJUGATE result (* Similar to a*i above *)ThuSj 
the correct NOTATION for this relation should be understood and written as 
c A | 

M’= <B | where the symbol "t" is the Greek symbol for DAGGER and is called the HERMITIAN CONJUGATION. 
Let us suppose that M | i > = | i' > and <j | M | i > = <j | i' > . If this relatation is TRUE, 
then <i|M f |j> = <i |j>> 

So the relation of two vectors < A | 

B > and < B | A > is once again a COMPLEX CONJUGATION or TRANSPOSITION . 

In our case the (i) and (j th ) matrix element of M f is the complex 

conjugate of the (j) and (i th ) matrix element of M. Thus, 

we can rewrite < i | M’ | j > = < j | M | i >* to indicate that there is an INTERCHANGE of (i) and (j) . 

A simpler relation is M f ij = Mij*. This is the notation from the Mathematician HERMITE, 

which is a REFLECTION about a DETERMINANT in the MATRIX 

and then proceeding with complex conjugation. For example, 

the Hermitian Conjugation of the matrix 

I 2 6+ M=( 2 ' 4_ ](* First reflect about the diagonal *) 

l4 -i 1 I \6+i 9 + i I 


Mi2 Mu 

M 2 2 M 23 


<*2 


a N 


= I 2 4+i 

\ 6 - i 9 + i 

(* Interchange the SIGN of the reflection and thus the complex conjugate as the final result *) 
The Hermitian Conjugate so far turns a BRA VECTOR into a KET VECTOR; and in terms of MATRICES, 
it flips (i) and (j); From the example above (6 + i) is TRANSPOSED to (4 - i), 
then complex conjugated to ( 4 + i) and (6 - i) respectively . 

The common notation ( Z = Z*) representing a REAL NUMBER 
is now replaced by HERMITIAN COMPLEX CONJUGATION OF MATRICES . So, 

M = M f and where the DIAGONAL ELEMENTS ARE ALWAYS REAL (4 and 6) 
while the COMPLEX NUMBERS ARE COMPLEX CONJUGATED ( + i and - l) 
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EIGEN VECTORS 


This can be defined as M | I > = | I > where (M) represents HERMITIAN and (I) represents EIGEN VALUES. 

If a matrix is Hermitian and the dimension of space is (n), 
then they will be n mutual ORTHOGONAL EIGEN VECTORS 
(* This defines a BASIS UNIT or BASIS VECTOR *) 

From the definition above, 

the EIGEN VALUES of HERMITIAN MATRICES are then REA1. So < I | M f = < I | X z * means 
An; is flipped and complex conjugated. So from the definition M | I > = An; | I 

< I | M = < I | Ai*. 

We can take the INNER PRODUCTS of both the equations and obtain, 

<I|M|I>=A I <I|I> and 

< I | M + | I > = < I | 

I>A;[* . If we substitute the similiraties of the two equations and cancel out <1 | I >, then 

Aj = Aj* = REAL. 

So , If a matrix is Hermitian, their Eigen Values will be equal to complex conjugates and they are REAL . 

Recall that we mentioned ORTHOGONALITY implies UNAMBIGOUSLY 
DISTINGUISHABLE. Let us now PROVE that different EIGEN VECTORS with DIFFERENT 
EIGEN VALUES are ORTHOGONAL, that is , they are NOT complex con jugates. 

Let M | I > = Aj | I > and 

M | I > = A] | I> . First, note 

< I | M = 3 | A] and so we have the two equations M | I > = A r | I > 

<3 | M =3 | A]. We can the multiply both the left_hand 
sides with (3) and (I) respectively to realize a substitution as we did earlier. So, 

< 3 | M | I > = A! < 3 | I > 

<3 | M f | I > = <3 | I > A 3 . So we have Aj = A : ; 

The requirement of the PROOF is that <3 | I > = 0. If we assume a value of 0 for the INNER PRODUCT, 
we obtain a BASIS VECTOR or an ORTHORNORMAL VECTOR. 

So far, we now have adequate MATHEMATICAL FORMALISM to INTERPRET the PHYSICAL MEANING OF 
THE STATE OF A SYSTEM. 

We begin with the notion that OBSERVABLES are represented by HERMITIAN OPERATORS where 
the POSSIBLE VALUES ARE EIGEN VALUES of (+1) and while EIGEN VECTORS CORRESPOND to the 
EIGEN STATES(Meaning that they are unambiguously distinguishable) 

The STATE of a SYSTEM is only OBSERVABLE when it is prepared or CONSTRUCTED using EIGEN 
VECTORS with EIGEN VALUES. 
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So , 

the PROBABILITY of X z , 

is represented by the SQUARE OF THE INNER PRODUCT of the VECTOR | A > with (I) . That is, 

Aj = < A | I > < I | A > (* The PROTECTION of |A> into the (I) DIRECTION *) 

So For ^ < A | I > < I | A> = 1 simply means that 

i 

< A | A > =1 which is the sums of the square of the components which is = 1; 
Thus the ctj = 1. 

i 

If we multiply the STATE VECTOR by a PHASE, for instance (e 1 ®) , 
then ALL PROBABILITIES ARE INVARIANT or UNCHANGED. This implies FEW independent COMPLEX NUMBERS. 


Let us now describe the ELEMENTS OF PAULI MATRICES which will be the 
EIGEN MATRIX VALUES applied as the OPERATOR representing the COORDINATES : 



So that an Eigen Matrix Value, 

multiplied by the Eigen Vector should = A Basis or Eigen Vector (* Orthonormal Vector*) That is, 

fon CT z : (l _*! ) (e ) = ( + {_ 0 i** 0 0 } } ) = ( 0 ) ' You can also show that CT * and are Eigen Matrices • 
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7. TIME EVOLUTION OF A SYSTEM 

CONSERVATION OF ENERGY AND SPIN 

We shall now analyze Mathematically, the notion of Linear Hermitian Observables. 


d I ifr > , 

The TIME DEPENDENT SCHRODINGER EQUATION is represented as i —= -i <H i/r > while 

at 1 

The TIME INDEPENDENT SCHRODINGER EQUATION is represented as <H | E ± > = E ± | E ± > 

(* The summation convention does not apply for the repeated indices of Ei *) Then, 

a | | 

— < i/r t L i/r t > = -l < i/r (LH - HL) 1 1/ > where 

at 1 1 

(L) represents AVERAGE VALUES 

(7Y) represents Hamiltonian 

(E) represents Eigen Vectors or Values 

a aL 

So — < L > = -i < [L, H] > where the shortened version is — = -i [L, H] 

at at 

meaning that the COMMUTATOR [L, H] of a REAL is ALWAYS IMAGINARY. 

We ' 11 now introduce the POISSON BRAKETS and the COMMUTATOR. 


The properties of Poisson Brakets is such that {A, B} = -{B, A}; So the COMMUTATOR [A, B] = - [B, A] . 

If we apply DIFFERENTIATION RULES for the PRODUCT {AB, C}, we obtain 
{AB, C> =A {B, C} + {A, C} B so the COMMUTATOR [AB, C] = A[BC] + [AC] Band where 
(AB, C - C, AB) = A (BC - CB) + (AC - CA) B 

thus the COMMUTATOR satisfies the ALGEBRAIC RULE of the Poisson Brakets. 

For the DIMENSIONAL STATES, {A, B} = - - [BA] can be rewritten as h {A, B> = -i [BA] 

h 

where h is the PLANCK 1 S CONSTANT. This implies the COMMUTATOR is negligibly small. 

So TIME DEPENDENCE is realized as | i/r t > = | j > where <Xj are the TIME_DEPENDENT COEFFICIENTS . 


Thus, 


dilr 

at 


y aa j(t) 

at 

3 


j > 


If we now insert the SCHRODINGER EQUATION to the equation above we write the DERIVATIVE as 

| *<t) > =X a 3< t > - i>H | =Z a 3( t ) - iw | £ i> 

3 j 

= ~ i 'Z J a Ut) * | Ej > 
j 

We know that the SUM of two VECTORS MUST have EQUAL COEFFICIENTS. So, 

a a,- 

-= -E-i a-i. This illustrates how the COEFFICIENTS change with TIME . 

at ] ] 

Therefore, the solution of the DIFFERENTIAL EQUATION of the COEFFICIENTS is an EXPONENT. That is, 
a j = a j (0) e lEj *• (* Recall that ^■= n [ x "~ 1 ] *) 


Now that we know the solution for the coefficients (a), 
we can then GENERALIZE the SCHRODINGER EQUATION at 

I <lft> = Zj a j<0) e " iE3t | j> 

where the coefficient (a) at time 
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(t 0 ) is the ENERGY DEPENDENT STATE (e _iEjt ) and is the PHASE. 

The PHASEj is the RELATION between ENERGIES and OSCILLATIONS. That is, the WAVE FUNCTION. 

So the sum of <i/r t | L | i/r t > =Zk a *k(0) e iEkt <k | L | j> Zj a j(0) e ‘ iEjt 

Note now from the GENERALIZED SCHRODINGER EQUATION that 

«At| = Zj i (0) e- tiE 3 t < j | 

. This is represented by the 

NOTATION (L 


Zjk “*k(0)“j(0) Lkje 1(Ek Ej)t is the HEISENBERG'S representation 
of the matrix formulation of the way AVERAGES CHANGE WITH TIME. 

This means that the two FREQUENCIES E k and Ej are associated 
with a MATRIX OBSERVABLE. 

So (L) CHANGES WITH TWO FREQUENCIES IN TIME. 
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Note that the MAGNITUDES of the COEFFICIENTS stay FIXED 

and the TIME EVOLUTION is associated ONLY with the RELATIVE PHASES. 

ID 

We shall now evaluate the relation n = — c 7 . 

2 

We can now apply the PAULI MATRICES mentioned to 
determine the STATE OF PREPARED SYSTEM in three coordinates (x, y, z) where 

id 

= -l[°x °zl ~ 

ID 

CT y = -i[a y CT z ] — and 

ID 

d z = -i[o z o z ] — = 0 (* You can work out the Matrices to prove its ORTHOGONAL ; 
Solution for working out matrices is^° ={ {ar+bt,as+bu},{cr+dt,cs+du}}*) 


The denominator (2) is arbitrary and will ease the substitution procedure later on . 
Therefor, the COMMUTATION RELATION of the SIGMAS (a) can be worked out as : 


/0 -i \ 

- 1 0 1 U 1 

i 0 

+1 1 

_ 2 I 0 1 1 

lie) 

l -l 0 I ' 

i -l 

0 ) 

2 l -1 0 ) 


2 i ct v 


Here we have taken the product of the commutators of the matrices. Recall 
to change the sign because of the order of the commutator. We thus get a + 
ve sum in the second step. (* Carefully work through the manipulation *) 

Recognize ( ^ g ) being related to a y which is j® j . If we multiply i by a y we obtain 
i ( ® ) = ( ® 1 g J . This NEATLY represents 2 ia y = [o z a x ] . Therefore we can ascertain that 

[° z cf x ] = 2 ia y 
[a x a y ] = 2 i a z and 
[a y a z ] = 2 ia x . 

We can substitute the values of the commutators to obtain : 
d x = -i[-2i] O y ^ = -CDO y 


. r . -j w 

d y = -i[2ict x J - 

and 

ct z = 0 


(DO x 


Its important to note that (i) disappers since (-i*-i=-l) and 
(-i*i = l) . This means that whem commutators have imaginary coefficients, 
the TIME DERIVATIVES definitively have a REAL outcome in general. 

Physically, This implies a PRECISION of ct x and a y along the ct z axis; 
similar to a gyroscopic precision or a rotor in a magnetic field in CLassical Mechanics. 

A 2D SPACE is the generalized HAMILTONIAN due to lack of independent operators. 


Note that u is the product of the magnetic field and magnetic moment of the spin. 
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0) 0) 

We can thus conclude that the ENERGY level of the system n-i = — a z is ALWAYS ± — . 

If it is flipped, we get w ENERGY. This means the SYSTEM is in a SUPERPOSITION STATE of either (±1) . 
If its - 1 there is NO RADIATION; while + 1 there ' sRADIATION. 

Thus, LINEAR OPERATORS act on EIGEN VECTORS and 
result in ORTHOGONAL EIGEN VECTORS where the AVERAGE VALUES 
must 0 for a SINGLE VECTOR STATE. 

For a COMPOSITE VECTOR STATE the system cannot be a seperate state and so the AVERAGE VALUE = 0 . 
This is a STATE of ENTANGLEMENT and is briefly outlined in the next chapter. 
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8. STATISTICAL CORRELATION AND FUNCTIONS OF VECTOR STATES 

What does it mean for an AVERAGE VALUE to be (-1) ? Here the situation is not of an inner 
product state. 

The state of a composite vector is mostly oppositely oriented where the Eigen 

vectors are CORRELATED statistically. Do not think of this in a Classical system 
otherwise youwill encounter INCONSITENCIES. This property cannot be simulated. 

The correlation of a composite vector occurs when the coefficients do not FACTORIZE . Reflect 
back to the construction of the commutator of <j z; So the product a z and z z = -1 
In a PRODUCT STATE there is NO CORRELATION. We can thus know information about a COMPOSITE SYSTEM but NOT 
the CONSTITUENTS. Its VECTOR STATE allows this nonsensible notion because it totally Mathematical. 

The low Energy of a state implies correlation . Sit down and penetrate 
through this by summing over a vector multiplied by a matrix. 

The DENSITY MATRIX for a PRODUCT or PURE STATE is =1. 

At what point does the interpretation become un_Mathematical ? Keep in mind the 
procedure is a DERIVATIVE of TIME . For instance, two connected computers cannot share 
information INSTANTENOUSLY. when you write a program to simulate a Vector State Space, 
accurate result simulation or the estimating precision is logically different. So 
dont be fooled. This means MEASUREMENTS are IRREVERSIBLE in UNITARY 
MATRICES. This is nature of simulators. Pretty intersting stuff so far. 


FUNCTIONS OF VECTOR STATES 


Let us now analyze OBSERVABLES, using PARTIAL DERIVATIVES. 
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Let us first define a FUNCTION where : 

i i • 50 

D 0 > —» 0 > So we can then say that D 0 (x) = —. Is it HERMITIAN ? 

1 1 dx 

The function can be written as an INTEGRAL . Observe. 

J d it r 5 0 i * 

0 * (x) — = 0 * (x) — 

dx J ax J 


The right_hand side can be interpreted 

r 50* 

as 0(x) -if it is Hermitian and is similar to the left hand side. 

J dx 


This is realized by INTEGRATION BY PARTS where the function or endpoints = 0. 

d ifr 

This is done by by SHIFTING THE DERIVATIVE — to 

dx 


3 0 * 
dx 


but we MUST now SHIFT THE SIGN of the left_hand side to get the result : 




50* 

0 (x) -. Notice that the equation is NOT EQUAL. 

dx 


There is a -ve, 

meaning the equation is opposite of each other. 


We can fix the situation by applying a ( -ve) Hermitian Operator such f-i — ] defined as ( P) . 

V dx ) 


So from equation (a) with (P) inserted we obtain : 

50 -j * 


(- 1 ) j*' <x) J' 


0*(x) 


dx J 


where (P) can represent MOMENTUM and written as ^ — J . Here h is to give the relation a sensible UNIT. 


Then the AVERAGE VALUE of < 0 | L | 0 > = 
dx) 




0* (x) 


0 (x) = REAL . The proof of this relation is in the integration by parts. 


This is illustrated by the shifting of the derivative and thus realizing a + ve (i) as shown below : 


-i JV (x) 


50 

dx 


(x) = +i 


‘J 


5 0* 
dx 


0 (X) . 


Thus j < 0 | L | 0 > = J 0 * (x) |-i ~ j 0 (x) = -i JV 


50 

(x) — (x) = +1 
dx 


‘J 


50* 

dx 


0 (X) . 


Observe the complex conjugation 

5 0* 

where 0* becomes - . You MUST CLEARLY see the complex conjugates. 

dx 

TWO VECTORS that are COMPLEX CONJUGATES of each other are REAL. 

We can now rewrite that : 

r 3 0 

< 0 I L I 0 > = -i 0 * (x) — (x) = REAL and thus is a HERMITIAN OPERATION and an OBSERVABLE. 

J dx 

The COMMUTATOR and the POISSON BRAKETS are related by the value (ifi) where 

[Xj P] = ill and {x, P} = 

5 

1 where (P) = -i — . Keep this relation in mind. It is the generalized Hermitian Operator. 

5 x 

For now, lets proceed with a Mathematical interlude. 
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FOURIER ANALYSIS 

Most FUNCTIONS can be represented in the form of a WAVE . That is, 

e lPx as Cosines and Sines where (P) is OSCILLATION and is modulated either as a short or long wavelength. 
Each wavelength has an AMPLITUDE which can be represented as i } (P) e lPx . So we can define the function 

i/r (x) = f aP i } (P) e 1Px . Note that a/2 tt is a CONVENTION. So, 

J V 2 7r 


i/r (P) = f i lr (x) e“ ipx . This is the BASIC FORM of the FOURIER THEORY where i/r (x) = i/r (P) , 

J * 


ax 
V 2 7r 

thus i/f (P) is the FOURIER TRANSFORM of i lr (P) with a (-ve) sign. 

This expresses the relation between (x) and (P) which is used to DEFINE POSITION and MOMENTUM. 


Let us now DEFINE the properties of the <5 FUNCTION. 

First consider the <5 function as an OPERATOR. Then we can write that : 

6 (x - x ') / (x) = /(x') <5 (x-x') since (x - x ') ^ 0. 

Thus the 6 (x - x ') = 1 illustrated by an AREA = 1. We can then represent the function where 


the 


J<5 (x-x 1 


) fW =/(*') =/(*') 


J 5 (x-x 


) dx. 


This is the 6 function in its entirety.Therefore, 


J 

J 


p co ap • . r r ax ~ 

(x ') = e lPx ilr (x) e“ lPx where we have inserted i/r (P) and simplified to obtain 

V 2 7 T J V2 7 T 

J rdP . 5P 

dx — e lP <* -x) ^ , Note that — e lP can be simplified by the notation A = 6 (x - x ') 

J 2 7 T 2 7 T 

So the FOURIER TRANSFORM is thus expressed as the SUM OF PLANE WAVES such that 


e -iP(x-x') =s (x _ x .) 


dP 
2 7T 

— e lx<p - p,) = <5 (P-P 1 ) ■ Obviously we have interchanged (P) with (x) . Importantly, 
2 7T 


note6 (x-x') and <5 (x ' - x) is a symmetrical function 
of its arguments such that the order of (x) and (x 1 ) does not matter. 


Now we can proceed and describe the notion of a particle on a line . 
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9. STATE OF A PARTICLE ON A LINE 


We can now describe in detail the meaning of CONSERVATION of Energy using the FOURIER 
TRANSFORM PAIR. 

Classically this is simply defined using (x, P) . In VECTOR STATE SPACE we must provide an ORTHOGONAL BASIS 
OF STATE (normalized or unit vector) so that we realize a LINEAR SUPERPOSITION of the VECTOR (BASIS) . 

In vector state space, the notation 

I *> 

| P > or 

| xP> represents too much information. Recall the formalism 

CTx 

CT y or 

We therefore cannot know both STATES simulatenously. This is 
called the HEISENBERG ’ s UNCERTAINTY PRINCIPLE. Its one or the other [x or P] . 

In Classical Physics a PARTICLE can be defined in both (x) and (P); 

While in vector state space the particle is described as one or the other . For instance; 

| ill > = J~3 (x) i/r (x) | x > means that the INNER PRODUCT of 

< x | 0 >= i/r (x) which is the projection of the WAVE FUNCTION and 

I It* (x) 1 A(X) =P(x) is the probability that (x) is at P (x) . (* VERY CLEAR AND NEAT *) 

Note that < x ' | x > =0 unless (x ' = x) . 

But now we are relating the density matrix for the INNER PRODUCT or PURE STATE. So, 

<x'|x>=l=<5(x-x')(* Recall that the Integral of the <5 function is = 1 *) 

So the TOATL PROBABILITY = 1. 

We have always described KETS as i It (x) and BRAS as i/r* (x) . That is | 
i/r (x) > and < 0* (x) | respectively. 

Thus, we can say , 

< 0 | 0 > = 6 (x - x ') and represented as 
i jr* (x) 0 (x) which when expanded is written as 

jr (x) 0 (x) dx 

If a particle is moving in 2 D, you cannot simulatenously know the x_coordinate of POSITION 
and the y_coordinate of MOMENTUM. This is the purpose of the FOURIER TRANSFORM PAIR. 

From the formalism of functions of vector states, we expect an EIGEN VALUE acting on an 
EIGEN STATE to give the same EIGEN VECTOR. So an OPERATOR (X) acting on a WAVE FUNCTION 
gives a number (X) times the wave function as shown below; 
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X 6 (x - x 0 ) = X 0 6 (x-x 0 ) where 6 (x - x 0 ) is the Eigen Vector and X 0 is the Eigen Value. 

This is similar to 

X i/r (x) = x i/f (x) and completes the formalism for (x) or POSITION. What about (P) or MOMENTUM? First, 

a 

note that (P) is an OBSERVABLE (REAL) and is CONSERVED. It is defined by the HERMITIAN OPERATOR -i —. 

ax 

. a 

So (P) = -ifi — . Thus, 
dx 

a i/r 

Pi/r (x) = -ih —. Now if we apply the Hermitian operation we obtain 
dx 

a <a p 0 

— = +i — i/r (x) . Recall that this is a partial derivative operation too. Then, 

ax h 

i/r (p) (x) = e 1 <x) means that the function of the wave 

is equal to the distance (x) has to move to change to 2 7r. In other words, 

— A = 2 7T or A = This a clear relation between MOMENTUM (P 0 ) and WAVELENGTH (A) . 

h P 0 


We can thus proceed and evaluate the INNER PRODUCT between two MOMENTUM STATES for | P D > . 
| P 0 > is represented by 

i/r Po (x) = e lf>0 <x) . It follows that the complex conjugation of the 


< P ' | P > 


T. 


,-iPo (x) e iP (x) 


ax 




X (P- P ' ) 


ex and thus factorized. Therefore, 


If P = P ' then the Integral of 1 becomes infinite. That is 1 dx = oo. The point is. 


In order to determine the NUMERICAL COEFFICIENT, 

we have to refer back to the 6 function and substitute to obtain 

!x(P- p ) = ( g (p _ p 1 ) anc j s i nce ih e the INNER PRODUCT of the 5 function is REAL, 




we have to divide e lxP by the V2 tx to obtain : 


<Ap„ (x) 


0 X XP o 


V2 7T 


. This is the WAVE FUNCTION of the particle with MOMENTUM (P) 


Thus, the WAVE FUNCTION of MOMENTUM can be illustrated 
< P | I/r > = I/r (P) whereas for POSITION is < x | i/r > = i/r (x) . (* i } 
the WAVE FUNCTION of POSITION can be illustrated as: 


as: 

implies an OBSERVABLE *) So 


e -i PX„ 

i/r Xo (P) = -and the solution is the 

V27T 



i/r (x) dx which is equivalent to the wave function of the FourierTransform. Thus, 


i/r (P) = i/r (x) which is exactly the wave function of FourierTransform. 

Therefore, the PROBABILITY of finding the wave function of momentum i/r (P) is : 
p (P) = i Y ( p) i/r (P) which is the square of the momentum 

(INNER PRODUCT) . This is similar to transposing from cs z to ct x Vector State.So, 

a 

If (P) =-i -—■ then 
dx 

. a 

(x) = + 1 — 

a p 

Note that TIME, TEMPERATURE, PRESSURE are not a Vector State Operators. 

An Electromagnetic Field and Angular Momentum are good examples of a CONJUGATE VARIABLE . 
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d i/r 5 i/r 

Any equation of the form — = -c — has the function 
at ax 

i/f (x - ct) . This solves the SCHRODINGER TIME DEPENDENT EQUATION . 

A particle with the COMPLEX FUNCTION i/r (x - ct) moves 
rigidly along the x_axis with a uniform velocity (C) ; 
where (c) represents the speed of light with the propertty of a + ve and - ve Energy when w = cP; 

Thus the PROBABILITY or AVERAGE VALUE or INNER PRODUCT of (x) moves with 
velocity (ct) . This is the interpretation of the particle in Vector State Space. 

Thus, the Hamiltonian representation of the motion of a particle can be written as : 

aw . J aw 

- =X and -= -P 

dP dx 

so that x = c and 

-P = 0 implying the particle is moving at the speed of light while MOMENTUM IS CONSERVED. 
(* This is similar to RECOIL, if something hits you. *) 


So far we can describe the HAMILTONIAN in terms of Vector State Space where : 

1 2 

T = Kinetic Energy = — mV 

P = Momentum = mV and 
P 2 

W = — = cP. 

2m 

a i/r 

Recall that Schrodinger ' s equation is represented as i — = -iw I i/r > 

at 

dii/ i . .a 

So ifi —= — (-i.fi) —i/r(x) . 
dt 2 m K ' dx 

(* P = -i-2- and (^) 2 =4; 

ex 'ax' ax 2 


the square of the 1st derivative is the 2nd derivative *) 

But (P) is squared ...So the equation then becomes: 

ifi^=^ (-ifi) 2 ^-i/r(x) which is described as a particle in one dimension moving 
along a line. Therefore, we can simplify the equation further to obtain: 
ifi^= — . This means that waves with different wavelengths move with 

dt 2m dx 2 

different velocities or different momenta move with different velocities. 


We conclude by noting that the basic rule for PARTIAL DERIVATIVES is that 
r F—dx = - r—G dx 

J dx j dx 

where we have shift the factor of the derivative at the cost of -ve sign. 
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10. UNCERTAINTY PRINCIPLE 


Let (Ax) 2 = J~i/r* i/f (x) x 2 ax = < i/r | x 2 | ^ > and 

2 f- 2 121 f a 2 i/r f , d dil/ 

(AP ) 2 = l/r(P) lA(P) P ap= < i/r P 2 Ur>= - dxi/r* (X) - = - axi/r* (X)- 

J 1 1 J ax 2 J ax ax 


where (P) = -i — thus (P = -- and so 

ax ' ' ax 2 


(AP)' 


'J 


d diJ/ 

— 1 1/* (x) —- Sx which is a complex conjugate and = REAL amd +vethus 
ax ax 


the INNER PRODUCT of (P) 2 . 

(* Recall the basic rules for Partial Derivatives and integrating by parts *) 


Consider the TRIANGULAR INEQUALITY (* the product of the magnitude 
of two vectors is bigger than the dot product between them and Cose t 1 *) 
where | A | 2 * | B | 2 > | A * B | 2 for the vectors : 

Substituting for (ax) 2 = | A | 2 and 

(AP) 2 = | B | 2 We would like to show that 

(AX) 2 * (AP) 2 >h 


We first asssume that Wave Function of (x) or i/r (x) is REAL 

. . d>!> 

| A > = X i/r (x) and | B > = Pi/r = -i —— 

ax 

If we considerthe INNER PRODUCT of the two vectors < A | B > 

and assume a REAL VALUE (since absolute values are REAL anyway, so we ignore the - ve signs), 
then 


<A | = | 


.2 


dij/ 


nx) — I ny) — <* noting *> 


1 par i par 

= -—(*)-— (y) 

2 J ax 2 J ay 

1 r dll/ 2 r dll/ 2 

-- ~ —— ( x ) — (y) • 

4 J ax J ax 


We can the INTEGRATE BY PARTS by shifting 
the derivatives keeping in mind we are expecting an Absolute Value. 
So, we have 


- f— a? r 

4 J ax J 


— ili 2 . (* noting — 

ay Sx 


1 = DERIVATIVE OF AN IDENTITY *) 


We know that the 


J* ! = 


which is the Integral of probability of all space if the wave function is normalized. Therefore, 

1 

each of the Integrals = 1 when the state vectors are UNITY. We are now left with the value —. 

4 

1 

Thus, (A x ) 2 * (A p ) 2 > — • Simplifying we have 
4 

(A x ) * (AP) > - and 

2 

a 

since (P) = -i — meaning (x) and (P) are NOT inverse to each other 

ax 
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( because their product cannot be a pure number ) — must be inverse to (x) . So , 

' ' h 

API h 

(Ax) * —— > - = (Ax) * (AP) > -• (* MIMIMUM UNCERTAINTY *) 
h 2 2 

This is the the Uncertainty Princple and is TRUE for all SPACE. 


We can now show that the TIME DERIVATIVE of MOMENT is equivalent to FORCE 

iTH I i/r > 


d a i/r 

where P = F = -i — . Consider the Schrodinger Wave Equation i — 

ax at 


Then, considering a particle moving in one Dimension along a line, we can show that 


. dil/(x t ) 


-4/ (X t ) +V (X) l/r (X t ) 


at 2 max 2 

which is the Schrodinger " s representation of how the Wave Function changes with TIME 
We can then get rid of the (i) on the left_hand 

side by multiplying both sides by (-i) and simplifying to obtain : 


diHx t ) 

at 

aiT (x t ) 


2 m ax 2 


i/r (x t ) - i V (x) i/r (x t ) since - i* i = 1 . Therefore, 


-i/r* (x t ) +iV (x) i/r* (x t ) which is a complex conjugate. 


at 2» ax 2 

We thus have both the Schrodinger 1 s Equations . 


Moving forward, the VELCOCITY OF A WAVE can be represented as : 

J r 

t P 

(X) * + J** (X) * 

: J'i/t* (x) i/r + J'i/t* (x) i/r which is a complex conjugate and (x) is REAL 


— < X > = — | i/r* (X) l/r (X) 

at a ■ 

-air 


So to get < x > 

we take twice the REAL part which is 2 (x) i/r since the imaginary part cancels. 

We then subtitute with the Schrodinger Equation 


ai jr(x t ) _ _i_ a 2 
at 2 m dx 2 

a 


i/r (x t ) - i V (x) i/r (x t ) to obtain : 


a r . r i a" i . . ... 

— < x > = 2 i/r (x)-i/r - l V (x) i/f noting that V (x) is imaginary .So 

at J L 2 m ax 2 - 1 


we cancel (-V (x) i/f) and the result 


= 2 fi/r* (X) -^*1 

J L 2 m ax 2 J 


We then must Integrate by Parts . The partial derivative operation yields : 

J r i d dilr-\ 

i/r* (x)-(* derivative with respect to x without changing time *) 

L 2 m ax ax J 

i r , r ai/ri i ei/f*ai/r 
2m J L ax J 2 m ax ax 

(*— = 1 and we must change signs *) The imaginary part cancels which then 

ax 


2— [V[— ]=-- [i/r* [— 1 • We know that -i —= P 
2 m J L ax ■> m J Lax J dx 


— = P = <i/r|P|</r> 


thus 


<<A 


i/r > which is the INNER PRODUCT of (P) divided by (m) . 


We have so far proved that 
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d < X > 

at 

where m 


: -and is the Classical relation between Velocity and Momentum 

m 

< x > . , 

-= < P > implies Mass * Velocity = Momentum. 

at 


That is one half of Classical Mechanics 


ap 


ax 
at' 


The othe half is — = F; which is the TIME DERIVATIVE OF MOMENTUM. 

at 


Note that 


ap 

at 


= P = F : 


av 

ax 


meaning that the time derivative of momentum is equivalent FORCE . 


a a r a , , 

— <P> = — U*(Xt) —I A (xt) -i 

at at J ax v ' 

ai/r 


—< i**—* 

ax 


J .* ei/r par 

i/r — + i -i/r* i/f. We can cancel the imaginary part so then 

ax J ax 

. rsiA* . ■ 

= + 1 — lA* lA 

J ax 

Observe that the INNER PRODUCT of < P > 

has a REAL and must be squared. We then substitute with the Schrodinger Equation to obtain : 


raiA* r 

i -\ 1 

a 2 x i 

— 

-1 — 

— -V (x) iA 

J ax L 

' '2m 

ax 2 J 


where 




1 a 2 x 

2 m ax 2 


J is imaginary and cancels . Thus, 


■; 


a i V_ 
ax 


V (X) !/r. 


V 2 m ax 2 


are now left with the equation 

-ai/r* 


can then be Integrated by parts and results in an imaginary number. We 

ai/r* 


the result 
thus 


par 

-V (x) i/r which must be complex conjugated as : 

J ax 

par r „ dt// 

= -V (x) i/r + i/r* V (x) -= REAL 

J ax J ax 

r dijf* 

V (x) -i/r + i/r* — which is the derivative-i/r 

J L ax axJ ax 

r r 3 lA* 1 

I V (x) ^ -iAj . If we now Integrate by parts, 

r * av . r . 

- hA (A — rewritten as - iA — iA 

J ax J ax 


av av 

= < lA* — lA > which is the INNER PRODUCT of — . 

ax ax 

We have so far proved that 

a av 

— < P > = - — in the sense of the AVERAGE VALUE in a smooth POTENTIAL in PHASE SPACE . 

at ax 

This means the WAVE FUNCTION is NOT disrupted or otherwise there is scattering . 


In conclusion, after the painful procedure but worthy. 

Below are the Wolfram Language CODE and INPUT to : 

(a) verify the calculation of MINIMUM UNCERTAINTY of a Radial Wave_Function of the Hydrogen Atom. 
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(* First Compute the energy eigenvalue from the differential equation *) 
(* The energy is independent of the orbital quantum number 1 *) 

'2r\i 2 
n 


4i[n_, 1_, r_] : = 


(n - 1 - 1) ! 


(n + 1) ! 


LaguerreL n - 1 - 1, 21 + 1 


■-] 
n 1 

. ,., r 1 I 2 i(i + i) 2 \ i 

energy[w_, r_, n_, 1_] := Simplify — D[w, r, r] t — D[w, r] - --- w + — w 

L wl r r 2 r;J 

Solve[ (energy[i/r[nj lj r], r, n, 1] - S//FullSimplify) == 0, 6] 

Out[a] := {{ S -> — JJ and 
It is REAL and + ve 


(b) Visualize Light Amplitude transmission through n layers of glass. 

Plot [ Evaluate [Table [l/ (l + k 2 ChebyshevU [n, Cos [k] + Sine [k] ] 2 ), {n, 0, 5} ] ], 
{kj 0j 10}, PlotRange-» All] 

Out[b] : 



(c) VISUALIZE and VERIFY the Quantum harmonic oscillator Wave Functions 
(* 

Input Wave Function 
Plot 

Normalize 

Compute the expectation value of x 2 
Verify that Momentum and Position wave 

functions for a Harmonic Oscillator have the same form *) 


iA [ n_, x_] : = 


HermiteH[n, x] e ~2 
V 2" n ! V^T 
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Plot [>[50, x] 
Abs[i/r[10, 


, {x, -10, 10}] 
x] ] 2 dlx 


x 2 Abs [i/r[10, x]] 2 dx 


{FourierTransform[i/r[2, x], x, o>], i 2 i/r[2, x] }//Simplify 
(* Notice that (i 2 ) is the OPERATOR*) 


Out[832] = 



Out[833] = 1 

21 

Out[834] = - 

2 

G ) 2 

e"T (l - 2 w 2 ) 
V2" 7T 1/4 


Out[835] = 


e-V (l-2x 2 ) , 

VT 7T 1/4 


(* Relate the output with 

— =- and —— where Position and Momentum are two frequencies in TIME *) 

at m at sx 
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